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Abstract 

'•^ ' We investigate the extreme value statistics connected with the dilute 

^ ^ ' Random Energy Model with integer couplings. New universality class is 

found. 

a' 

Extreme value statistics [1] statistics play an important role in statistical me- 
chanics and theory of probability [2-5]. They describe the minimal and maximal 

j^ , bounds for the distribution of a large number of independent random quanti- 

ties, and it is the main reason of their recent applications in statistical physics 
of disordered systems. 
' O ■ Usually one considers sequence of AI = 2^ independent random variables Ei 

S ! with the same distribution 

^ ■ piE) =< 6{E, -E)> (1) 

Cn ■ Let us find minimal and maximal members Emin and Emax of our sequence, then 



> 

o 



consider their distribution at the limit of large N. When p{E) is independent 
from N, at TV — > oo limit only 3 classes of universalities has been found [1] for 
the distribution of p(-E'rnin). After rescaling of a; = E^i^ one has for them: 

O' Pi(a;) = l,a;>0,Pi(a;) = l-exp[ 

a\'- P2(a;) =0,a;<0,P2(a;) = l-exp[ 
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C^ ■ Pi{x) = l-exp[- 

y ' Now we are going to consider another (physically interesting) case, when our 

p{E) is a superposition of N random variables. 
fH ' In the present paper we will investigate extreme value statistics which are gen- 

Q , erated by the dilute-coupling REM. Surprisingly, we got new types of them due 

O ■ to physical peculiarities of the model. 

The dilute coupling REM is described by the hamiltonian 

>< ■ ^ 

b. ^ = ^2-^ J»i..*pSi,..s,^ (5) 

(ii..ip) = l 

where only aN quenched bonds are nonzero, and they take values ±1 with equal 

probability. Si = ±1 arc the Ising spins. 

In the limit of large p there is a factorization in probability of different energy 



levels and we have 2^ independent random variables (energies) with the same 
probability distribution p{E). For our choice of integer couplings we have for 
distribution p{E). 

p(E) = -^ / dfce^^('=)-'=^ (6) 

where 

f{k) = Q:lncosh(fc) (7) 

This representation is correct for 

1^1 < a (8) 

Otherwise p{E) is strictly zero. 

For the case a < 1 we should keep in mind, that there are 2"^ different classes 

of energy configurations with 

equal energies in each one. Energies from different classes are distributed inde- 
pendently according to (6). 

At low temperatures the system is frozen in minimal energy configuration (chozen 
among the all 2^ ones). Thats why the physics of REM is connected closely 
with extreme value statistics problem. 

As known, several phases are possible in the considered model. For a > \ there 
is a phase transition to the spin glass phase and at high temperatures the sys- 
tem is in the paramagnetic phase. For the case a < 1 only the paramagnetic 
phase is possible. The limiting case a ^ 1 corresponds to the border between 2 
phases. 

We hope, that different universality classes are possible, connected with those 
phases. Let us consider 2^ independent random variables with (6), where ^p^k) 
is even analytical function with respect to k. It is just enough for the most of 
situations. For the case of hamiltonian (5) one has, that 

lim ip{k) = a{k - ln2) + o(exp(-2fc)) (10) 

k — *oo 

We are considering only distributions , which have this asymptotic form, or 
when corrections are polynomial 

lim Lp{k) ^ a{k - In 2) + o{k-^) (11) 

k — >oo 

We are going to consider now the extreme value statistics for 2^ random vari- 
ables with our representation for p{k). 

Let us define the probability 0{E), as the integral of these expressions in the 
interval [—oo,E]. We have 

-t f-ioo 

9{E) = / dfce"'^W-'=-^//c (12) 



Here our integration loop overpasses the point from the left size. 
Let us find the distribution of the minimal bound. We take one variable Ei in 
our interval [—00, E], the others to right. As the first one could be chosen from 
the all 2^, we have an expression [4] 

PiE) = 2^p{E){l e{E)f-^ = -^(1 - ^(^))'" (13) 

We need to study the expression 9{E) in the region, when it is small. Here 

(1 - 6l(S))2" « exph2^6'(£;)] (14) 

Its density is nonzero only near the E = E^- We need to find asymptotic of 
function 9{E) near the E^, defined by condition 

2^0(i?c) = 1 (15) 

For the different choices of a we have to consider different asymptotic regims 
for 0(E). In this way different universality classes arose. 
While calculating (12) for the saddle point we have an equation 

'P'(fc) = -| (16) 

There are three possibilities for Ec- The first one 

\E,\<aN (17) 

has the same statistics, as in the case when (p{k) goes to infinity faster than lin- 
early (this is the case of ordinary REM [4] with normal distribution of energies). 
Two others are possible also: 

\E,\ ~ aN (18) 

\E,\ > aN (19) 

In the last case we should take complex solutions of (16). Let us first calculate 
the first case (18). Two expressions for the function 0{E) coincide (at the limit 
of large N) and we have 

9(E) = 1 ^Nv(k)-kE (20) 

ky/2TTN(p"{k) 

One can write the expression for the desired probability distribution for the 
minimum of our 2^ random variables: 



k^2TrN^"{k) " J J^ Nip(k)-kE /2-|^^ 



dE ky^ 271 Nip"{k) 



For our asymptotic form (15) transforms into 

2N 



ky^2TrNip"{k) 



^Nv(k)-kE, = (22) 



where k is defined from (16). For the E < Ec our probability is almost one, and 
above it is almost zero. For a small deviation 

u= {E- Ec)kc < N (23) 

we obtain Gumbcl distribution 

P{u) = cxp(?i — cxp(u)), (24) 

where u changes in the interval — oo, oo. Let us now consider the statistics for 
limiting situation, when a ^ 1. We shall consider the distribution (7), but the 
corresponding limiting statistics is the same for the whole class of exponential 
corrections. At the saddle point 



tanh(/c) = E/{aN) (25) 



we have 



200 



where 



e{E) = —l e"^lncosh(fe)-fcAr/^_ (26) 

^ COsh(fc) ^^aNg(E/(aN) 

1 + X 1 + X 1 — X 1 — X , ^ 

5(^)-^-ln^- + ^-ln^- (27) 

Let us do notation 

N,=Na,^ = l-e (28) 

We are considering the case, where -^ — 1 = ei ^ 0. In the limit e ^ we have 

OiE) = ^-^^ exp (-iV, m 2 + ^ ln(2/.)) (30) 

Let us check, when such a behaviour is correct. We should take condition 



Now let us consider the region Ne ~ 1. Of course, at the saddle point we cannot 
approximate our integral by the gaussian one for our values of e, since the second 
derivative tends to zero. We should directly calculate integral. 
For the complex A: = fci + ik2 we have 

Incosh(fc) = coshfci + l/21n[l - (1 - tanh^ fci) sin^ ^2] 

+i arctan[tan fc2 — (1 — tanh/ci) tanfc2] (32) 

This is an exact expression. Let us consider the case /ci ^ 00. In this limit 



In cosh(A:) = fci - In 2 + e^"^^^ , tanh fci = 1 - 2e 



-2fei 



So we have to choose the position of fci , while performing integration 



1 



ki^ioo 



^a7V[lncosh(A;)-(l-e]fc 



/k 



(33) 



(34) 



-ln(^) 
2 ^e' 



The saddle point equation gives the choice 

fci 
With such a choice of fci we get 

iVi Incosh(fc) - Ni{ki + ik2) + Nie{ki + ifc2) 

wiV(e/2 + e/21n 2N / e - e sin fc2 + ie (fc2 - sin fc2 cos ^2 ) ) 

Accuracy of this expression is exp(— 2fci) ^ e. Our expressions are valid in a 
region 



(35) 



(36) 



(37) 



We get 



e{E) wcxp( -7Viln2+— ln(- 



F^2\ — 1 „ — e sin k'. 



2 fci cos(e(fc2 



— / dk2{kl + ^2) 

TT 

— sinfc2Cosfc2)) + fc2sin(e(fc2 — sinfc2COsfc2)) 



For the limit of small e last expression transforms into 

Ne 
e{E) =cxp(-]Viln2 + — ln(2/e)) 

We have for the desired distribution for u = (^ — l)iV 



P(u) 



- 1 + in — I exp 
u 



uln ■ 



N 



exp u In 



N 



(38) 



(39) 



(40) 



Such distribution is correct in the interval — oo < u < 0. What about nega- 
tive values of u? For the choice (7 )with integer aN we have, that there is a 
strict bound for energy distribution equal to TVa, so for that case we can put 
zero probability for negative values of u . This expression should exist for any 
bounded discrete distribution of energies. 

In [5] authors consider fully connected version of REM, here perhaps could be 
only classic distribution among the (2)- (4). The case of bounded continuous 
distribution brings to (2)- (3). 
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